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Amendments to the Claims: 

This listing of claims will replace all prior versions, and listings, of claims in the application: 
Listing of Claims: 

1 . (Currently Amended) A ligand screening apparatus which screens for a ligand 
that binds to a protein when coordinate data of the protein of single chain or plural chains is 
given, the apparatus comprising: 

a post-structural-change protein coordinate data selecting unit that e ff e cts 
conducts structural change in consid e ration of in the coordinate data of the protein while 
considering dynamic behavior , wherein said structural change is performed using an induced- 
fit parameter reflecting induced fit on the coordinate data of the protein and s e l e cts post- 
structural -change protein coordinate dat a is selected ; 

a spatial point designating unit that designates a spatial point at which 
superposition with the ligand is to be conducted from the post-structural-change protein 
coordinate data selected by the post-structural-change protein coordinate data selecting unit; 

an interaction function calculating unit that calculates an interaction function 
when the protein and the ligand bind to each other using the spatial point designated by the 
spatial point designating unit and a ligand coordinate data of the ligand; and 

a ligand evaluating unit that evaluates the ligand that binds to the protein based 
on the interaction function calculated by the interaction function calculating unit. 

2. (Previously Presented) The ligand screening apparatus according to claim 1, 
wherein the interaction function calculating unit calculates the interaction function using 
Sscore(i j) shown in Formula 1 : 
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A. 

Sscore(i,j)=Y, 



when i is not equal to j 
«*[exp {- fe; }-/?}/ 

7k +d f) : 

when i is equal to j 



[Formula 1] 



wherein dj/ is a distance between i-th spatial point and j-th spatial point in the target 
protein, djj C is an interatomic distance between i-th atom and j-th atom in a compound, a is a 
coefficient for making Sscore(i j) the maximum value when the group of spatial points in the 
target protein and the compound completely overlap with each other, p is a coefficient for 
giving a threshold value by which it can be defined as "overlapping". 

3. (Previously Presented) The ligand screening apparatus according to claim 1, 
wherein the interaction function calculating unit further comprises an interaction function 
optimizing unit that carries out optimization so as to make the score of interaction function 
maximum. 

4. (Previously Presented) The ligand screening apparatus according to claim 3, 
wherein the interaction function calculating unit further comprises: 

an interaction energy optimizing unit that calculates interaction energy with 
the protein for the superposed ligand after optimization of the interaction function by the 
interaction function optimizing unit, and optimizes the interaction energy while finely 
adjusting conformation of ligand 3D structure data. 

5. (Previously Presented) The ligand screening apparatus according to claim 4, 
wherein the ligand evaluating unit further comprises: 

a reevaluating unit that causes execution of the interaction function calculating 
unit after largely changing conformation of ligand 3D structure data following optimization 
by the interaction energy optimizing unit, and reevaluates the ligand that binds to the protein 
based on the interaction function calculated by the interaction function calculating unit. 
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6. (Previously Presented) The ligand screening apparatus according to claim 1, 
wherein in calculation of any one of the induced-fit parameter and the post-structural-change 
protein coordinate data or both, the post-structural-change protein coordinate data selecting 
unit calculates a normal mode for the protein coordinate data, determines intensity of 
fluctuation of each amino acid, and conducts a molecular dynamic calculation using the 
intensity of fluctuation as a constraint condition. 

7. (Currently Amended) The ligand screening apparatus according to claim 6, 
wherein the post-structural-change protein coordinate data selecting unit calculates a 
fluctuation value of a dihedral angle of a main chain according to normal mode calculation, 
and conducts a molecular dynamic calculation while setting the fluctuation value as a 
coefficient of force K in the molecular dynamic calculation shown by Formula 2 or Formula 
3: 



wherein Erot represents energy of a dihedral angle of a main chain atom in a 3D 
structure of a protein, § represents a dihedral angle of the main chain atom, <|)0 represents a 
standard value of the dihedral angle of the main chain atom, when a value of Krot (a 
coefficient of force) is large, (j> is constrained by (|)0[[.]] a 



wherein Epos represents position energy of the main chain atom in a 3D structure of a 
protein, r represents a coordinate of the main chain atom, rO represents a standard value of the 
coordinate of the main chain atom, when a value of Kpos (a coefficient of force) is large, r is 
constrained by rO. 

8. (Withdrawn) The ligand screening apparatus according to claim 1, wherein the 
interaction function calculating unit uses the interaction function to which a dynamic property 
function representing dynamic property of protein is added as "elastic energy". 




[Formula 2] 




[Formula 3] 
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9. (Withdrawn) The ligand screening apparatus according to claim 8, wherein the 
interaction function calculating unit adapts "U collision" as "elastic energy" which is a 
function shown by Formula 4 in consideration of local flexibility of the protein 



M N 

U collision =Z!X^( / >>) 

(p (i, j) = Kcollision * (Rcollision (/, j)-R) 2 

[Formula 4] 

wherein M is a number of atoms in an active site that prohibit collision, N is a number 
of atoms of the ligand, interatomic distance R between an i-th atom of a main chain or a side 
chain with a little dynamic behavior in an active site, and a j-th atom in the ligand is not more 
than collision distance "Rcollision(ij)", <J)(i j) is calculated. 

10. (Previously Presented) The ligand screening apparatus according to claim 1, 
wherein the interaction function calculating unit uses the interaction function to which a 
normal mode analysis result or secondary structure determination result of the protein is 
added as a dynamic property function that represents a dynamic property of the protein. 

1 1 . (Currently Amended) A ligand screening method which screens for a ligand that 
binds to a protein when coordinate data of the protein of single chain or plural chains is given, 
the method comprising: 

a post-structural-change protein coordinate data selecting step that e ff e cts 
conducts , with a computer, structural change in consideration of in the coordinate data of the 
protein while considering dynamic behavior , wherein said structural change is performed 
using an induced-fit parameter reflecting induced fit on the coordinate data of the protein and 
s e l e cts post-structural-change protein coordinate data is selected ; 
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a spatial point designating step that designates a spatial point at which 
superposition with the ligand is to be conducted from the post-structural-change protein 
coordinate data selected by the post-structural-change protein coordinate data selecting step; 

an interaction function calculating step that calculates an interaction function 
when the protein and the ligand bind to each other using the spatial point designated by the 
spatial point designating step and a ligand coordinate data of the ligand; and 

a ligand evaluating step that evaluates the ligand that binds to the protein 
based on the interaction function calculated by the interaction function calculating step; 

wherein the interaction function calculating step calculates the interaction function 
using Sscore(i ,j) shown in Formula 1 : 

x 

Sscore(i,j)=Y, 

ij 

wherein dy S is a distance between i-th spatial point and j-th spatial point in the target 
protein, dij C is an interatomic distance between i-th atom and j-th atom in a compound, a is a 
coefficient for making Sscore(i J) the maximum value when the group of spatial points in the 
target protein and the compound completely overlap with each other, p is a coefficient for 
giving a threshold value by which it can be defined as "overlapping". 

12. (Cancelled) 

13. (Previously Presented) The ligand screening method according to claim 11, 
wherein the interaction function calculating step further comprises an interaction function 
optimizing step that carries out optimization so as to make the score of interaction function 
maximum. 

-6- 

WASH_7649359.1 



when i is not equal to j 
ax[exp {-fa-d<f}-p\ 
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14. (Previously Presented) The ligand screening method according to claim 13, 
wherein the interaction function calculating step further comprises: 

an interaction energy optimizing step that calculates interaction energy with 
the protein for the superposed ligand after optimization of the interaction function by the 
interaction function optimizing step, and optimizes the interaction energy while finely 
adjusting conformation of ligand 3D structure data. 

15. (Previously Presented) The ligand screening method according to claim 14, 
wherein the ligand evaluating step further comprises: 

a reevaluating step that causes execution of the interaction function calculating 
step after largely changing conformation of ligand 3D structure data following optimization 
by the interaction energy optimizing step, and reevaluates the ligand that binds to the protein 
based on the interaction function calculated by the interaction function calculating step. 

1 6. (Previously Presented) The ligand screening method according to claim 1 1 , 
wherein in calculation of any one of the induced-fit parameter and the post-structural-change 
protein coordinate data or both, the post-structural-change protein coordinate data selecting 
step calculates a normal mode for the protein coordinate data, determines intensity of 
fluctuation of each amino acid, and conducts a molecular dynamic calculation using the 
intensity of fluctuation as a constraint condition. 

17. (Currently Amended) The ligand screening method according to claim 16, 
wherein the post-structural-change protein coordinate data selecting step calculates a 
fluctuation value of a dihedral angle of a main chain according to normal mode calculation, 
and conducts a molecular dynamic calculation while setting the fluctuation value as a 
coefficient of force K in the molecular dynamic calculation shown by Formula 2 or Formula 
3: 

Erot = Krot (<|> - <|>0) 2 [Formula 2] 



-7- 

WASH_7649359.1 



Atty. Dkt. No. 079088-0102 



wherein Erot represents energy of a dihedral angle of a main chain atom in a 3D 
structure of a protein, § represents a dihedral angle of the main chain atom, §0 represents a 
standard value of the dihedral angle of the main chain atom, when a value of Krot (a 
coefficient of force) is large, <j> is constrained by (^[[.Jk 



wherein Epos represents position energy of the main chain atom in a 3D structure of a 
protein r represents a coordinate of the main chain atom, rO represents a standard value of the 
coordinate of the main chain atom, when a value of Kpos (a coefficient of force) is large, r is 
constrained by rO. 

18. (Withdrawn) The ligand screening method according to claim 11, wherein the 
interaction function calculating step uses the interaction function to which a dynamic property 
function representing dynamic property of protein is added as "elastic energy". 

19. (Withdrawn) The ligand screening method according to claim 18, wherein the 
interaction function calculating step adapts "U collision" as "elastic energy" which is a 
function shown by Formula 4 in consideration of local flexibility of protein: 



<p (i, j) = Kcollision * (Rcollision (/, j) - R) 2 

wherein M is a number of atoms in an active site that prohibit collision, N is a number 
of atoms of ligand, interatomic distance R between an i-th atom of a main chain or a side 
chain with a little dynamic behavior in an active site, and a j-th atom in the ligand is not more 
than collision distance "Rcollision(ij)", <{)(i,j) is calculated. 

20. (Previously Presented) The ligand screening method according to claim 11, 
wherein the interaction function calculating step uses the interaction function to which a 
normal mode analysis result or secondary structure determination result of the protein is 
added as a dynamic property function that represents a dynamic property of the protein. 




[Formula 3] 



M N 



U collision =SX^( / >^) 



[Formula 4] 
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21 . (Currently Amended) A non-transitory computer-readable medium having a 
program which makes a computer execute a ligand screening method which screens for a 
ligand that binds to a protein when coordinate data of the protein of single chain or plural 
chains is given, the method comprising: 

a post-structural-change protein coordinate data selecting step that e ff e cts 
conducts structural change in consid e ration of in the coordinate data of the protein while 
considering dynamic behavio r, wherein said structural change is performed using an induced- 
fit parameter reflecting induced fit on the coordinate data of the protein and s e l e cts post- 
structural-change protein coordinate dat a is selected ; 

a spatial point designating step that designates a spatial point at which 
superposition with the ligand is to be conducted, from the post-structural-change protein 
coordinate data selected by the post-structural-change protein coordinate data selecting step; 

an interaction function calculating step that calculates an interaction function 
when the protein and the ligand bind to each other using the spatial point designated by the 
spatial point designating step and a ligand coordinate data of the ligand; and 

a ligand evaluating step that evaluates the ligand that binds to the protein 
based on the interaction function calculated by the interaction function calculating step; 

wherein the interaction function calculating step calculates the interaction function 
using Sscore(i j) shown in Formula 1 



wherein dy S is a distance between i-th spatial point and j-th spatial point in the target 
protein, dy is an interatomic distance between i-th atom and j-th atom in a compound, a is a 



when i is not equal to j 




[Formula 1] 



when i is equal to j 
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coefficient for making Sscore(i j) the maximum value when the group of spatial points in the 
target protein and the compound completely overlap with each other, P is a coefficient for 
giving a threshold value by which it can be defined as "overlapping". 

22. (Cancelled) 

23. (Previously Presented) The non-transitory computer-readable medium having a 
program according to claim 21, wherein the interaction function calculating step further 
comprises an interaction function optimizing step that carries out optimization so as to make 
the score of interaction function maximum. 

24. (Previously Presented) The non-transitory computer-readable medium having a 
program according to claim 23, wherein the interaction function calculating step further 
comprises: 

an interaction energy optimizing step that calculates interaction energy with 
the protein for the superposed ligand after optimization of the interaction function by the 
interaction function optimizing step, and optimizes the interaction energy while finely 
adjusting conformation of ligand 3D structure data. 

25. (Previously Presented) The non-transitory computer-readable medium having a 
program according to claim 24, wherein the ligand evaluating step further comprises: 

a reevaluating step that causes execution of the interaction function calculating 
step after largely changing conformation of ligand 3D structure data following optimization 
by the interaction energy optimizing step, and reevaluates the ligand that binds to the protein 
based on the interaction function calculated by the interaction function calculating step. 

26. (Previously Presented) The non-transitory computer-readable medium having a 
program according to claim 21, wherein in calculation of any one of the induced-fit parameter 
and the post-structural-change protein coordinate data or both, the post-structural-change 
protein coordinate data selecting step calculates a normal mode for the protein coordinate 
data, determines intensity of fluctuation of each amino acid, and conducts a molecular 
dynamic calculation using the intensity of fluctuation as a constraint condition. 
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27. (Currently Amended) The non-transitory computer-readable medium having a 
program according to claim 26, wherein the post-structural-change protein coordinate data 
selecting step calculates a fluctuation value of a dihedral angle of a main chain according to 
the normal mode calculation, and conducts the molecular dynamic calculation while setting 
the fluctuation value as a coefficient of force K in the molecular dynamic calculation shown 
by Formula 2 or Formula 3: 



wherein Erot represents energy of a dihedral angle of a main chain atom in a 3D 
structure of a protein, <|> represents the dihedral angle of the main chain atom, §Q represents a 
standard value of the dihedral angle of the main chain atom, when a value of Krot (a 
coefficient of force) is large, (J) is constrained by (|)0[[.]] i 



wherein Epos represents a position energy of the main chain atom in 3D structure of a 
protein, r represents a coordinate of the main chain atom, rO represents a standard value of the 
coordinate of the main chain atom, when a value of Kpos (a coefficient of force) is large, r is 
constrained by rO. 

28. (Withdrawn) The program according to claim 21, wherein the interaction 
function calculating step uses the interaction function to which a dynamic property function 
representing dynamic property of protein is added as "elastic energy". 

29. (Withdrawn) The non-transitory computer-readable medium having a program 
according to claim 28, wherein the interaction function calculating step adapts "U collision" 
as "elastic energy" which is a function shown by Formula 4 in consideration of local 
flexibility of protein: 




[Formula 2] 




[Formula 3] 



M N 



U collision = ('">./) 



<p (i, j ) = Kcollision * (Rcollision (z, j)-R) 2 
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[Formula 4] 

wherein M is a number of atoms in an active site that prohibit collision, N is a number 
of atoms of ligand, interatomic distance R between an i-th atom of a main chain or a side 
chain with a little dynamic behavior in active site, and a j-th atom in the ligand is not more 
than collision distance "Rcollision(i ,j)", (|>(i j) is calculated. 

30. (Previously Presented) The non-transitory computer-readable medium having a 
program according to claim 21, wherein the interaction function calculating step uses the 
interaction function to which a normal mode analysis result or secondary structure 
determination result of the protein is added as a dynamic property function that represents a 
dynamic property of the protein. 

31. (Cancelled) 
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